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ABSTRACT 

Halos, filaments, sheets and voids in the cosmic web can be defined in terms of the 
eigenvalues of the smoothed shear tensor and a threshold Ath- Using analytic methods, 
we construct mean maps centered on these types of structures for Gaussian random 
fields corresponding to cosmological initial conditions. Each map also requires a choice 
of shear at the origin; we consider three possibilities. We find characteristic sizes, 
shapes and other properties of the central objects in these mean maps and explore how 
these properties change with varying the threshold and smoothing scale, i.e. varying 
the separation of the cosmic web into different kinds of components. The mean maps 
become increasingly complex as the threshold Ath decreases to zero. We also describe 
scatter around these mean maps, subtleties which can arise in their construction, and 
some comparisons between halos in the maps and collapsed halos at final times. 


1 INTRODUCTION 

The web of cosmological large scale structure is both strik¬ 
ingly evident in observations and simulations and difficult 
to characterize. Its existence has been known for several 
decades (Zel’dovich, Einasto & Shandarin 1982; Shandarin 
& Zel’dovich 1983; Einasto et al. 1984; Bond, Kofman & 
Pogosyan 1996), and a huge number of different approaches 
have been used to study it since. Recent overviews of many 
of these directions can be found in the Proceedings of the 
Tallinn conference (2014). 

These cosmological structures evolved via gravitational 
instability from initial conditions which have been character¬ 
ized by measurements of the cosmic microwave background 
(e.g., Planck Collaboration (2013)). These early time ini¬ 
tial conditions are well described by Gaussian random fields, 
which can also be classified into cosmic web structures. An 
assortment of different approximate methods can be used 
to estimate properties of the late time counterparts of the 
structures found in the initial conditions. Properties include 
the fraction of mass which will collapse into a certain type of 
cosmological structure, the number of such structures, etc. 
Classic works include Press & Schechter (1974); Bardeen, 
Bond, Kaiser & Szalay (1986); some textbook descriptions 
are found in Peacock (1998) and Mo, van den Bosch & 
White (2010). For instance, peaks in the initial conditions 
can be associated with collapsed halos at late times on aver¬ 
age (Kaiser 1984), and properties of the peaks can be con¬ 
nected to final time properties of the halos (for instance by 
tracing back particles in simulations to their early time con¬ 
figurations such as in, e.g., Dalai, White, Bond & Shirokov 
(2008); Robertson, et al (2009); Ludlow and Porciani (2011); 


Elia, Ludlow & Porciani (2012); Despali, Tormen & Sheth 
(2013); Borzyszkowski, Ludlow & Porciani (2014); Ludlow, 
Borzyszkowski & Porciani (2014)). Halos are the most stud¬ 
ied cosmological structures, in part due to their close iden¬ 
tification to galaxies (Blumenthal, Faber, Primack & Rees 
1984). 

This note extends studies of cosmic structures within 
the specific cosmic web classification considered in Alonso, 
Eardley & Peacock (2014). (Their detailed particular defini¬ 
tions are in §2 below.) In that work, they compared analytic 
predictions for initial random Gaussian fields (equivalently, 
their values linearly extrapolated to final time) and mea¬ 
surements of fully evolved dark matter simulations at late 
times (redshift 2 = 0.1). For some values of their parame¬ 
ters, including the primary one which we focus upon in this 
work, the volume fractions of the structures for Gaussian 
random fields and structures for fully evolved dark matter 
simulations counterparts at late times were comparable. The 
analytic mass fractions followed the same trends as found in 
the final time simulations, but were somewhat too high for 
sheets and too low for filaments. (Halos, which have the 
largest overdensities and thus might be expected to have 
the largest differences in the two methods, contribute a very 
small fraction to both the mass and volume fractions in this 
structure classification.) 

This analytic classification of cosmological Gaussian ini¬ 
tial conditions has further consequences, as the power spec¬ 
trum also determines mean spatial properties around a given 
central point. We analyze these spatial properties in this 
note. We find mean properties of regions around different 
choices of origin, classified as halos, filaments, or sheets, re¬ 
gions we call mean halos, mean filaments and mean sheets. 
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(The boundary around a “mean void” is more challenging to 
define as for many parameter choices the mean background 
is also a void, so that a boundary is not evident). The cen¬ 
tral shear values of these mean objects are an additional 
parameter which must be chosen. We explore three natural 
choices for these mean object central shear values, and mea¬ 
sure sizes, shapes, shear profiles, and mean linearly extrapo¬ 
lated overdensities of the resulting mean halos, filaments and 
sheets. We also discuss scatter around these mean configu¬ 
rations, subtleties in their construction, and some possible 
relations between mean halos and final time collapsed halos. 
We consider structures in dark matter only. 

Background is given in §2: the particular shear-based 
web classification we use, analytic estimates for average 
shear in different objects, the construction of mean maps 
and the measurement of some of their properties. In §3 we 
consider three choices of central shear, showing mean map 
examples for different objects and trends in their proper¬ 
ties when changing the parameters defining the cosmic web. 
We also discuss scatter around the mean. In §4 we turn to 
possible connections of the mean halos to collapsed halos 
at late times. In §5 we summarize and conclude. The three 
appendices give details of the calculations (shear, Zel’dovich 
displacement and scatter), subtleties which arise in some in¬ 
terpretations of one of our construction methods and prop¬ 
erties of mean maps corresponding to slightly different con¬ 
straints at the origin used by Papai and Sheth (2013). 
The power spectrum used is for a ACDM universe with 
{nm,^bC,h,(78,n) = (0.31,0.022,0.68,0.83,0.96) (Planck 
Collaboration 2013), calculated with the Eisenstein & Hu 
(1999) transfer function. 


2 DEFINITIONS AND METHODS 

The cosmic web can be separated into halos (sometimes 
called nodes), filaments, sheets and voids. We will use the 
term halo and node interchangeably, as both terms are com¬ 
monly used for these structures. There are many different 
web classifications based upon density, morphology, shear 
and beyond, some of which concentrate on just hnding one 
particular type of structure (e.g. halos or voids) and oth¬ 
ers which classify an entire map. We use a classification of 
the latter sort, assigning a type of cosmic structure to every 
region based upon its shear. 

2.1 Shear 

We begin with details of our shear definition and classifi¬ 
cation. We follow Desjacques and Smith (2008); Papai and 
Sheth (2013) for dehnitions and for much of the analysis 
methodology below. The shear is defined on the initial Gaus¬ 
sian perturbations as 

Cy(r) = ^didj^ir) . (1) 

Our potential 'l?(r) has the same spatial dependence as the 
Newtonian potential (j> at early times but is time indepen¬ 
dent: 4>(r) = a)/A tvGC pbD{a). Here G is Newton’s con¬ 

stant, pi, is the mean background density of the universe, and 
D{a) is the linear growth factor describing how linear struc¬ 
ture grows as the universe expands with scale factor a. The 


derivatives di are with respect to the comoving coordinates 
r. The potential ‘I>(r) is related to the initial density per¬ 
turbations evolved linearly to the present time, (5(r), via the 
Poisson equation: 

^af<I.(r) = 5(r) . (2) 

i 

At times earlier than the present time the linearly evolved 
overdensity obeys 

5L{r,a) = D{a)5{r) . (3) 

The linearly evolved overdensity at the present time, (5(r) at 
position r is defined in terms of (the linearly evolved to the 
present time) density via p(r)/p — 1. 

The factor a is the square root of the integral of the 
power spectrum P{k) associated with 3{r), smoothed on 
scale Rs 

(7^ ^ ^ J dkCPik)WhAk) , ( 4 ) 

i.e., the root mean square fluctuations of the density contrast 
S{r). We use a Gaussian window function, 

WRAk) = e-i^^^^^\ (5) 

in part to connect more closely with Alonso, Eardley & 
Peacock (2014), as we use their parameters as a start¬ 
ing point below. One of the parameters of the classifica¬ 
tion is the smoothing scale Rs- We will show results from 
Rs = 4:h~^Mpc, studied in depth by Alonso, Eardley & 
Peacock (2014), but note trends we saw for other cases we 
explored, Rs = (0.5,1, 2, 8.9,10,19) h~^Mpc. 

Again, to make the notation clear, because D{a = 1) = 
1, 5{r) equals the linearly evolved density perturbation to¬ 
day and a corresponds to its commonly used dehnition as 
the square root of the integral of the power spectrum linearly 
evolved to the present day. Some other approaches include 
D{a) in the overdensity, and thus in 4>. As the fluctuations 
a are associated with the overdensity, in our definition they 
have a factor of D{a) as well, and so ^ij remains time inde¬ 
pendent. This is one reason for introducing a into our def¬ 
inition of shear. Other definitions of shear also sometimes 
take out the trace, again not the case here. When the shear 
in numerical simulations is considered, Eq. 1 is often used as 
well, but with the the fully evolved nonlinear potential (de¬ 
termined by the evolved densities in the simulation), rather 
than its linearly evolved counterpart used in the analytic 
work in this note. 

The potential 4? can also be used to Hnd the Zel’dovich 
approximation to the displacement of initial positions at a 
later time. In the Lagrangian approach to structure forma¬ 
tion, the final position x of an initial mass component with 
initial position q is related via the displacement 4'(q, t): 

x = q+^(q,f). (6) 

In the Zel’dovich approximation, 

= D(a)V,-I> . (7) 

For the initial conhguration, with coordinate q, multiplying 
the linear hnal time density by D{a) and taking a small 
gives the initial density (which is linear because it is small, 
5i(r, a) — D{a)S{r) ~ 0 as a —>■ 0 and x —>■ q). 

The shear can thus be thought of either as a property of 
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the linearly evolved final density perturbations, or as a prop¬ 
erty of positions at early times (which then change, using 
Eq. 6 ). We will discuss properties in both interpretations. 


2.2 Cosmic web classification 

The 3x3 shear matrix ^ij has three eigenvalues, traditionally 
ordered as Ai > A 2 > As. Their sum, the trace of the shear, 
is proportional to the density by Eq. 2, 

tr Cij = (8) 

^^ <7 

i 

Halos, filaments, sheets and voids can then be defined as 
regions where (3,2,1,0) respectively of the Ai are above some 
threshold Ath- In one way of estimating later structure for¬ 
mation, an initial region is expected to collapse at late times 
if it linearly evolves to cross some overdensity threshold Sc- 
As the linearly evolved overdensity 5{r) = Xi{r), and 
regions classified here as halos might be expected to collapse 
at late times (this is an assumption), this relation suggests 
some minimum Ai, i.e. nonzero Ath, for initial conditions. 
Note that this condition also means that the mean back¬ 
ground, Ai = 0, is also classified as a void. There are thus 
two scales in this classification, the smoothing Rs and the 
threshold Ath. Many different choices have been used in the 
literature. 

The same eigenvalue classification is sometimes used for 
dark matter simulations, but (as mentioned above) using 
the shear associated with the fully evolved nonlinear po¬ 
tential, see, e.g., Hahn et al (2007); Forero-Romero et al 
(2009) and references therein. For example, in dark matter 
nonlinear simulations evolved to redshift zero, Hahn et al 
(2007) found Rs = 2h~^Mpc and Ath = 0 separates struc¬ 
ture into these different classes in a way that matches vi¬ 
sual estimates, while Forero-Romero et al (2009) considered 
Rs = 0.88h~^Mpc and Ath = 0 , 0 . 2 ,0.4., !., 2 . and found im¬ 
provements for Ath > 0 based upon volume fractions and 
mass fractions. 

One comparison between the above analytic and simu¬ 
lation approaches is by Alonso, Eardley & Peacock (2014). 
In both approaches, they measured the mass fraction and 
volume fraction in the four cosmic structures as a func¬ 
tion of Ath, and for two values of Rs- They found reason¬ 
able correspondence of volume and mass fractions for fila¬ 
ments, sheets and voids between the analytically calculated 
linearly evolved perturbations and the numerically evolved 
fully nonlinear simulations. Given this, we looked for further 
characteristic properties of the cosmic web in the analytic 
construction and report these properties here. We focus for 
the most part on the parameters used in Alonso, Eardley & 
Peacock (2014).^ 


^ There are other structure classifications using the eigenvalues, 
for instance based upon their sum, 5(r) by Eq. 2, where each 
type of structure corresponds to (5(r) crossing some particular 
threshold 5th (Shen, Abel, Mo & Sheth 2006). The distributions 
of halos using Ath vs 5th was done in Lam, Sheth & Desjacques 
(2009) where it was found that the Ath constraint worked well. 
In analytic work inspiring the work here, Papai and Sheth (2013) 
considered a mixed constraint, using both Ath and 5th to define 
halos and voids, and compared density anisotropies at large dis¬ 
tances between the analytic (linearly evolved to final time) and 


We will use these definitions of halos, filaments, sheets 
and voids to analytically calculate mean configurations of 
each, for initial random Gaussian field configurations with 
different choices of Rs and Ath. Unlike Press-Schechter the¬ 
ory where one changes the smoothing scale to consider dif¬ 
ferent mass halos, here the one smoothing scale produces a 
separation of a full configuration into halos, filaments, sheets 
and voids, i.e. one smoothing scale completely classifies an 
entire map. Our mean objects correspond to one such spe¬ 
cific smoothing scale and threshold. 


2.3 Average shear eigenvalues 


The analytic definition of halos, filaments, sheets and voids 
in terms of Ath allows the calculation of volume and mass 
fractions in each type of object and the average value of 
each of the eigenvalues Ai,A 2 ,A 3 in the Gaussian random 
field initial conditions. 

The probability distribution for the ordered eigenvalues 
Al > A 2 > A 3 of a Gaussian random field is (Doroshkevich 
1970) 

J3(Ai, A 2 , As) = —e ^ (Al — A 2 )(Ai — A3 )(A2 — As) 

SvS-TT 

( 9 ) 

where /i = Ai -|- A 2 -t- A 3 and I 2 = A 1 A 2 -t- A 2 A 3 -i- Ai A 3 . Note 
that our A; differ by a factor of a from that in Doroshkevich 
(1970) because we include a in our definition of shear, Eq. 1. 
Expressions for pairs of eigenvalues or individual eigenvalues 
are found in Lee and Shandarin (1998) and another inter¬ 
esting rewriting is in Alonso, Eardley & Peacock (2014). 

The probability of the eigenvalues restricted by some 
condition C, given for example in terms of Ath, is given by 

P{C) = f d®Ap(Ai,A 2 ,A 3 ) . (10) 

■ic 

One case of interest for us is where the shear in the region 
C obeys the constraints for a halo, filament, sheet or void 
in terms of Ath. The mean eigenvalues for this particular 
constraint C, which we will denote by Xi^cA, are 

>xp(Ai,A2,A3)-p^ . 

( 11 ) 

In this example, for each structure the corresponding nor¬ 
malization P{C), Eq. 10, has the same integration limits 
as the the corresponding integral over p(Ai, A 2 , A 3 )Ai (it 
changes for different structures). For example, for a halo. 


halo 

filament dAi fff dX 2 dXs 

sheet rfAi/_"-dA 2 /_"LrfA 3 

f"ZdX^ffldX2ffldX, 


\ halo 

'^i,CA 


/ a °1 f^tl dX3p{Xi,X2X3)Xi 

/aTh '^A 3 P(Ai, A2A3) 


( 12 ) 


simulated (nonlinearly evolved to final time) calculations. The 
mixed constraint makes the identification of filaments and sheets 
less direct. We give properties using the Papai and Sheth (2013) 
constraint in Appendix §C. Structures classified using eigenval¬ 
ues of the Hessian of the density (e.g. Pogosyan et al (2009)) and 
the Hessian of the velocities (e.g. Hoffman et al (2012)) have also 
been considered in great detail by many authors. 
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We use the notation \i,cA to indicate these constrained av¬ 
erages. It depends both upon Ath and the type of structure 
(halo, filament, sheet or void). In practice, for any given Ath 
we integrated these expressions numerically using mathe- 
matica. 

Another constraint goes beyond restricting eigenvalues 
in terms Ath for each structure, requiring the <1> from which 
they come to also correspond to extremal values of the po¬ 
tential = 0). In this case the shear obeys conditions to 
be a halo, filament, sheet or void in terms of Ath, and its po¬ 
tential also obeys di^ = 0. The measure for the volume av¬ 
erage of points obeying this constraint is found in Bardeen, 
Bond, Kaiser & Szalay (1986), starting with the probability 
distribution for the random Gaussian fields <&, di^, ^ij times 
the constraint on di$. They were interested in peaks and 
minima in particular, as their main interest was for halos 
and voids (they also had a restriction on $, which we do not 
impose here, and thus it drops out of the expressions). Their 
initial measure for integrating the probability density over 
volume, shear and di^ becomes an integral over the eigen¬ 
values similar to that in Eq. 11, except with an additional 
measure factor IA 1 A 2 A 3 I (from changing variables in the con¬ 
straint Soidi^O) from di^ to position and integrating over 
position). 

/ dV / d®C / d^d^5{r - rp)p{^ij, . 

= /d3A/d3a<I><5(9$)p(Ai,A2,A3)|AiA2A3| . ^ ’ 

Here the integral again is over the values of Ai obeying the 
constraint to be a halo, hlament, sheet or void. We will call 
the corresponding average shear values Ai,extr. So, for exam¬ 
ple, for a halo, one has: 

^halo ^ -Ch A 2 A 3 )|AlA 2 A 3 |Ai 

/aTh -f^th '^•^3P(Ai, A2A3)|AiA2A3| 

(14) 

These two kinds of average shear, \i,cA and Ai,extr, are 
shown in Fig. 1 below for the 4 kinds of objects, as the 
dashed and solid lines, and are quite similar. 

We will consider one other shear average below, that 
of the mean objects we construct. The eigenvalues of the 
average shear within these mean objects will be called Ai. 

2.4 Mean shear maps: defining objects and 
characteristic properties 

One can go beyond the analytic volume fractions, mass frac¬ 
tions and the constrained shear eigenvalue averages for each 
structure type by using the power spectrum. The power 
spectrum can be used to calculate mean spatial properties 
of the shear around any given central shear configuration. 
We make maps of these mean shear values around a given 
central shear value, and study properties of regions within 
them, in particular those regions sharing the same cosmic 
web classification as the central point. 

We use the method of Papai and Sheth (2013), follow¬ 
ing Bardeen, Bond, Kaiser & Szalay (1986), to find the mean 
values of ^ij(r) (and thus Ai(i')) everywhere around a given 
configuration (0) — ^ij at the origin. Their calculation, 
for our case of interest, is summarized in Appendix §A. De¬ 
noting the constraint at the origin as Co , we thus get a map 
(iij ( r)|C'o), a mean shear at every point. This shear map de¬ 
pends not only on Go as explicitly indicated, but also upon 


the linear power spectrum and smoothing. The shear eigen¬ 
values at the origin corresponding to the origin constraint 
Go are denoted as A? 

Once the mean shear (Cij(r)|Go) is known everywhere 
around a given central configuration, a choice of Ath sepa¬ 
rates the mean map into halos, filaments, sheets or voids. If 
the origin constraint Go corresponds to a halo, the nearest 
points are usually also halo points, similarly for filaments, 
sheets and voids. We consider the central objects in the mean 
maps in the most detail. The central object in a mean map 
is the region sharing the same cosmic web classification (e.g. 
hlament) as the origin and which are continuously connected 
to the origin via regions with this same classihcation. Far 
away from the origin, the correlations and thus the mean 
shear go to zero and the mean conhguration becomes the 
mean background (again, this is a void region for Ath > 0). In 
between the origin and mean background, depending upon 
parameters, several different types of structures may occur. 
For instance mean halos tend to have hlaments and then 
sheets around them before reaching the mean background.^ 
Similarly, hlaments tend to have sheets around them (and 
sometimes halos within, towards their ends). Voids do not 
always have a clear boundary because for Ath > 0 they also 
can include the mean background. 

Several quantities can be associated with these result¬ 
ing mean map objects, such as volume, shape, and average 
values of shear XA. All of these depend upon Go, Ath and 
Rs. 

We report the volume of these mean objects in terms of 
the Lagrangian volume for a halo associated with the same 
window function and smoothing scale Rs, Vkaio Lag. This La¬ 
grangian volume is found by by integrating the background 
average density at early times (as SL{r,a) — >■ 0 ) against 
the window function to give M{Rs) and then converting to 
14aio Lag = M{Rs)/pb- For our Gaussian window function, 
Vhaio Lag = (2tvA^'^RA The correspondence between shear 
and initial Gaussian perturbations means that one can view 
the volumes of objects in our maps as Lagrangian volumes 
as well, as their volume also directly gives their mass at early 
times. 

For one choice of shear at the origin. A® = Ai,extr, 
boundary conditions for the Zel’dovich displacement at the 
origin are also specified, as 9;$ = 0. In this case we also cal¬ 
culate the mean linearly extrapolated Lagrangian displace¬ 
ment at late times given by Eq. 7. The mean displacement 
of the edges of the mean objects gives one way of estimating 
(using linear evolution) their final size. 

An additional measurement we make is of the aver- 


^ For simplicity we will always take a coordinate system such 
that the shear tensor is diagonalized at the origin, with the largest 
value of shear along x, second along y and smallest along 2 . 

3 See the nice analytic toy model in Alonso, Eardley & Peacock 
(2014) for an intuitive example of structure around a halo. In 
terms of the shear eigenvalues, configurations going from all three 
eigenvalues above threshold to no eigenvalues above threshold 
would generally be expected to pass through configurations with 
two or one of the eigenvalues above threshold. 

The average eigenvalues are calculated by taking all points clas¬ 
sified as belonging to the structure by their eigenvalues, averaging 
their shear over the volume of the object, and then finding the 
eigenvalues of this average shear. 
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lengths), average shear eigenvalues within, Xi, linearly ex¬ 
trapolated average overdensity(5), and in some cases mean 
Zel’dovich linearly extrapolated displacement . Most of 
our focus is on the objects which include the origin (objects 
not connected to the origin will be discussed briefly in Ap¬ 
pendix §B). 

We show examples below of cross sections of these re¬ 
gions along the 3 pairs of spatial axes and of the three shear 
eigenvalues and linearly extrapolated overdensity along the 
3 spatial axes. 


3 SHEAR MAP EXAMPLES 

The mean maps for a given cosmology are determined by the 
choice of central shear {^ij |Co) = (Ai,A 2 ,A 3 ), the thresh¬ 
old defining the separation into different structures, Ath, the 
smoothing Ra, the window function and the power spec¬ 
trum. We hold the window function (Gaussian) and power 
spectrum fixed and show results for Rg = Ah~^Mpc (com¬ 
menting on changes for other values). The threshold Ath and 
central shears Co vary as indicated. 


Figure 1. Average values of 3 constrained eigenvalues as a func¬ 
tion of threshold Ath for each kind of structure. The solid line 
is Xi CA and the dashed line is Ai^extn the dotted line is the 
threshold Ath- Eigenvalues for halos (red), filaments (blue), sheets 
(cyan), and voids (yellow, with black central lines for clarity) are 
shown at upper left, upper right, lower left, and lower right. These 
curves are also the central values A? for two of the three mean map 
constructions in §3.1, A? = Xi cj^ for choice 2 and A? = At^extr 
for choice 3. The filled squares are for choice 1, where A? is fixed 
to make the average shear Xi within the central object satisfies 
Xi = Xi CA - The mean maps hereon use the same colors (for voids, 
yellow alone) for each kind of structure. 

age of the linearly extrapolated Eulerian overdensity at fi¬ 
nal time in these mean objects, given by (S) — Xi. In 
simulations, this has been measured by running the clock 
backwards, starting with collapsed halo regions in nonlin¬ 
ear simulations. Particles in the final halo are tagged and 
then traced back in time, and the average of the shears of 
these particles in the initial conditions is used to calculate 
the linearly extrapolated Eulerian overdensity at final times. 
Despali, Tormen & Sheth (2013) found for lower final halo 
masses that on average dextrap ~ (5c(l+0.2cr(J?s)), and Dalai, 
White, Bond & Shirokov (2008) found similar trends for 
higher mass halos.® We quote overdensities in terms of this 
reference 5extrap below. 

To summarize, for central halo, filament and sheet re¬ 
gions in mean maps we will measure volume, shape (axes 

® The definition of halo at final times and the cosmological pa¬ 
rameters differed in these studies, and the measured relations had 
a lot of scatter, so these numbers are only approximate, another 
similar estimate was found in Elia, Ludlow & Porciani (2012). See 
§4 for more discussion of relations between collapsed halos and 
halos in this shear based classification. 


3.1 Choices for shear at origin, A°, of mean maps 

To construct a mean map, a choice of shear at the origin, Co, 
is required. Mean maps can be constructed for any choice 
of shear at the origin Co, but for a given type of structure 
and Ath, we found three choices of particular interest, all 
motivated by the analytic calculations of average shear given 
some constraint. Two are associated with Ai,cA, found via 
Eq. 11, and one with Ai,extr, which includes Eq. 13 in the 
calculation of the average shear eigenvalues. They all yield 
a value for A^, the shear at the origin as follows (implicitly 
for choice 1, explicitly for choices 2,3). 

1 . Xi = Xi^cA : 

Origin shear choice 1 fixes A® implicitly, by requiring that 
the average Xi within a particular mean object in the 
map, Xi, corresponds to Xi^cA for the corresponding type 
of structure. That is, the mean object with shear A° at 
the origin has average shear corresponding to the average 
analytic shear, Eq. 11 for that kind of structure. For 
example, for a halo in the mean map, A(*“*° = AJ'cai where 
the left hand side is the average over the values in the mean 
map, and similarly for the corresponding filament (using 
Af,cA) etc.). In practice, we change the central value of 
shear until the mean shears within the halo region centered 
at the origin (which changes in size and shape as the central 
shear changes), obeys Eq. 12. One possible interpretation 
is that this gives an “average” such object (e.g. halo), in a 
smoothed simulation. Practically, we could not always find 
solutions A° satisfying the constraint Xi = Xi,cA for all Ath 
and all structures, in addition, other complications in this 
implicit definition can arise (Appendix §B). Our examples 
below are in the regime where the associated corrections 
from these possible complications seem small enough to be 
ignored and where we could find solutions. 

2. Xi = Xi^cA ■ 

Origin shear choice 2 sets the central shear A^ equal to 
the average shear Xi^cA in Eq. 11 (similar to Papal and 
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Figure 2. Cross sections of a central mean halo (top left), filament (top right) and sheet (bottom) along three pairs of axes. For each 
central object, A? is calculated for the chosen Ath (listed in each plot) using choice 1: the average eigenvalues Ai = \i^cA for the 
corresponding type of structure (halo, filament or sheet). At top left, the mean halo (red) is surrounded by a filament region (dark blue), 
then a sheet region (cyan) and finally a void (yellow, note mean background is also a “void”). The smoothing Rs = Ah~^Mpc for all 
three, the threshold Ath is as indicated. The largest eigenvalue is along x. For the central object in each, volume and linearly extrapolated 
overdensity (5) are also listed. )No mean void is illustrated, as shear corresponding to a void at the origin can smoothly transform to 
the mean background, which is also a void. It is thus difficult to define the volume of the void and find a central value using this choice, 
choice 1, for A?.) 


Sheth (2013), who used it for long distance properties 
at final times and whose work inspired this study). This 
means putting the average shear for a certain structure 
(for instance a halo) at the origin and then measuring 
properties of the object made up of the enclosing region 
of points sharing the same cosmic web classification. 
In more detail, for a halo, the origin would have shear 
eigenvalues A® = ^^c°a- Although slightly more complicated 
to interpret, these maps are much easier to construct, as 
the central shear is explicitly and immediately calculable. 
Mean properties around such a point might be interpreted 
as mean properties around an average point in the corre¬ 
sponding structures in a large map, rather than properties 
around the central point of a structure (i.e. the peak of a 


halo does not have average halo properties).® 

q \0 _ \ 

/'i — '^i,extr ' 

Origin shear choice 3 is the volume average of shears obeying 
some particular cosmic web constraint where in addition il? is 
at an extrema, i.e. where fttl? = 0 for all i. This is reasonable 
to expect for a center of a halo or void (minimum or peak) 
and we can also consider this condition for the center of a 


° The mean objects constructed this way might have sizes more 
closely related to, e.g., how far (on average) one might have to 
go, from an average point for a given structure, before the object 
ends. The approximate radial position of such an average point in 
the mean maps above for A^ = Xi^cA can be read off from Fig. 6. 
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Figure 3. Cross sections analogous Fig. 2, but for central shear A? = Ai^extr (case 3 in §3.1). Colors, Rs, and Ajj, as in Fig. 2. In addition, 
initial edge positions displaced using the mean Zel’dovich approximation, Eq. 7, are connected across the origin via black lines (if they 
do not go beyond the origin). No mean void is shown as the central configuration smoothly transitions to the mean background, which 
is also a void. 


filament or sheet. These A® = Agextr are calculated with the 
extra measure factor derived from Eq. 13. Intuitively, the 
mean objects corresponding to these central shear values, 
at least for halos and voids, might correspond to stacking 
objects in a map, centered on their extremal points (which 
is often the point used to define the center of such objects). 


In Fig. 1 we show the central values for halos, filaments, 
sheets and voids for these three constructions, as a func¬ 
tion of Ath- As central value 1 required implicit solving of 
an equation, only a few values are shown, as points (and 
again, we could not find solutions for all Ath in this case, 
and the construction is ill-defined for voids). The solid line 
is for case 2, A; = \i,CA, the integral in Eq. 11, the dashed 
line is for case 3, the volume average over extremal points, 
A° = Xi^extr- The central values are very similar for the three 
constructions. 


3.2 Mean object cross sections 

With these central values we can calculate the mean shear 
around the origin, and classify each region as halo, filament, 
sheet, or void. In Fig. 2, we show example cross sections 
of regions obeying constraint 1 for the central object, i.e. 
where Ai within the central object is approximately equal to 
the corresponding AgcA, along the three pairs of axes. We 
use colors to distinguish between structures: red for halos, 
dark blue for filaments, cyan for sheets and yellow for voids.^ 
The mean background (Ai = 0) is also a void (yellow) for any 
Ath > 0. (As voids comprise the eventual mean background 

^ For the central halo and filament, we use parameters studied 
in detail in Alonso, Eardley &: Peacock (2014), smoothing Rs = 
4:h~^Mpc, and threshold Ath = 0.762. They chose Ath by trying 
to get similar analytic volumes in each of the different types of 
structures. For the sheet shown we changed Ath to 1.25 to avoid 
subtleties discussed in Appendix §B. 


© 0000 RAS, MNRAS 000, 000-000 






















8 Aung and Cohn 


into which these structures are found, and the majority of 
the volume, as noted earlier it was not clear how to impose 
Ai = \i,CA for their case. Thus no void example is shown.) 
As anticipated, the mean central halo (at top left. Fig. 2) is 
surrounded by a hlament, then a sheet, and then eventually 
reaches mean density (which is a void), similarly the mean 
hlament (at top right. Fig. 2) is surrounded by a sheet and 
then void and mean density, and so on. The mean hlament 
does not extend between halos for this value of Ath although 
as Ath is lowered halos do eventually appear in the mean 
hlament endpoints. The mean halo has comparable size in 
all three directions, the hlament is extended in one direction 
compared to the other two, and the sheet has one direction 
much smaller than the other two, as expected. The larger 
directions are a few times larger than the smoothing scale 
Rs. 

Fig. 2 also lists the mean extrapolated overdensity and 
volume for each central object. (By construction, for choice 
of origin 1, (<5) = a Xi,CA-) For the Ath and Rs here (0.762 
and 4h~^Mpc respectively), the halo volume is slightly 
smaller than expected for the Lagrangian region associated 
with a halo of mass M{Rs), and the extrapolated density is 
higher, V ~ O.GVjialo Lag, ^ ~ 1.45extrap. 

The change in these quantities using a different choice 
of Ai can be seen in Fig. 3, where we show mean halo, fil¬ 
ament and sheet maps for the same Rs and Ath values, but 
choice 3 for shear at the origin, A^ = Ai,extr. They are simi¬ 
lar to their counterparts in Fig. 2, but slightly smaller, as is 
expected with their smaller values of A° (see Fig 1). Choice 
3 of shear at the origin also Hxes 9$ = 0 at the origin, 
implying values for the mean Zel’dovich displacement for 
points in these structures. In Fig. 3, black lines along each 
axis connect the position of the initial edges, after they are 
shifted by their mean Zel’dovich displacements (V$). (No 
line is shown if the displacement crosses the origin. A dis¬ 
placement crossing the origin means that the estimates of 
the Hnal edge values are well outside of the range of their 
validity.) Under these displacements, the mean halo con¬ 
tracts and the direction with the largest shear at early times 
becomes the most contracted at late times, as expected and 
seen as well in Despali, Tormen & Sheth (2013), see also Lud¬ 
low, Borzyszkowski & Porciani (2014); Robles, Dominguez- 
Tenreiro, Onorbe & Martinez-Serrano (2015). In addition, 
the mean filament and sheet both contract (only slightly in 
some directions), while the void, not shown, expands in all 
directions. As mentioned earlier, the final time estimate for 
the linear displacement of the initial region edges provides 
one (very rough) estimate of the mean region’s final size in 
the Lagrangian picture. ® 

More generally, the shapes of the mean objects in these 
maps as the threshold Ath changes can be seen in Fig. 4, 
for Rs = 4:h~^Mpc and our 3 choices of shear at the origin. 



0.5 1 


Figure 4. The extent of halos, filaments and sheets, along the 
x,y,z axes when A? = A^ (heavy solid lines) and A^' = At extr 
(light dashed lines), as a function of Ath- Filled squares are coun¬ 
terparts for halos, filaments and sheets, at the thresholds shown, 
when Ai = Xi^cA, case 1. The smoothing scale Rs = Ah~^Mpc. 
The largest eigenvalue is along x. Examples of objects for case 1 
and case 3 are shown in Fig. 2 and Fig. 3 respectively. 


(Note that as these are the maximum of each coordinate 
that the actual filament length is, for instance twice the 
maximum z value.) Solid points give the extent along the 3 
axes for values of Ath where we solved choice 1 (At = Xi^cA), 
solid lines denote choice 2, mean maps centered on A°Ai,cA, 
and the lighter dashed lines are for choice 3, A® = Ai,extr. 
We only sampled central values for condition 1 because it 
required solving an implicit equation for each Ath and each 
different type of structure, and again, we could not always 
find a solution. 

For a fixed smoothing scale (here Rs), raising the 
threshold Ath in general raises the average shears within ob¬ 
jects, Xi. As can be seen, this simultaneously decreases the 
volume of the associated mean region for any object beside a 
void. The shapes of halos, i.e. ratios of extent along different 
axes, change strongly as the threshold Ath varies.^'' 


® This doesn’t always happen for other parameters. For instance, 
for the other set of Rs, Ath explored in detail by Alonso, Eardley & 
Peacock (2014), Rs = 10h~^Mpc, Ath = 0.261, the trend reverses, 

U ^ 2Vhalo Lag, ^ ^ 0-5^extrap- 

® A final overdensity estimate can also be made from the ratio 
of the final to initial volume in this picture. This does not always 
work well in our examples as the edges of structures be decreased 
beyond their original distance from the origin (implying the ap¬ 
proximation has broken down long before this point). 


For other smoothings Rs, not shown here, we found that the 
sizes of the structures along the different axes (in units of Rs and 
relative to those seen in Fig. 4) increase for lower Rs < Mpc 
and decrease for larger Rs > Mpc. For example, at low Ath, 
the size along a given coordinate axis {x, y, z) for a given Ath is 
about twice as large (2/3 as large), in units of Rs, for Rs = 0.5 
(19.0) h~^Mpc. The ratio of largest sheet edge to largest filament 
edge is close to constant as Rs changes, for Ath > 0, but the 
filament/halo largest edge ratio changes (drops) with increasing 
Rs until Ath 0.2. 
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For choice 3 of shear at the origin, the trends of the 
mean Zel’dovich displacements of the edges of the mean ob¬ 
jects could also be studied. We found that the halo edges 
contract for all Ath and that as Ath increases the displace¬ 
ments go through the origin and become negative, indicating 
the approximation is well beyond its range of validity. For 
hlaments and sheets, although 2 or 1 of the edges are/is 
displaced quickly towards the origin even for small Ath, for 
the other edge or edges, as Ath increases, the displacement 
either is small at first (although causing the displaced edge 
to move towards the origin), or causes the object to first ex¬ 
pand (the two larger sheet axes) before again having the dis¬ 
placed edge move towards the origin. The mean Zel’dovich 
displaced position of the edges divided by the original posi¬ 
tion of the edges takes roughly the same form for all choices 
of smoothing we studied. 

For many values of Ath and Rs, the mean objects around 
average points resemble those in Figs. 2,3 above, albeit with 
different side lengths as just noted. However, as Ath tends 
to zero, not only do central structures get much larger as 
seen in Fig. 4, but also qualitatively different behavior ap¬ 
pears. Cross sections for the case Ath = 10“'* and case 2 , 
A° = A i,CA, are shown in Fig. 5. The mean maps are shown 
for a much larger range of scales, because unlike the two 
earlier cases shown in Figs. 2,3, where outside the central 
regions the shear quickly tends to a void, for this low Ath 
non-void structures persist out to large distances. (We con¬ 
sider Ath = 10“* rather than exactly zero as the functions 
used to calculate the shears have limited accuracy.) For this 
threshold, in part because voids and halos only comprise 8 % 
of the total volume each, with sheets and hlaments taking 
42% each, the structures which are not the mean background 
(not voids) are extremely large. In addition, unlike Fig. 2, 
mean hlaments have halos within their endpoints, mean 
sheets have hlaments within their outskirts, mean voids have 
sheets outside of them, and so on, as can be seen. (The pres¬ 
ence of halos within the ends of hlaments is characteristic 
of individual hlaments in simulations, which are often found 
by looking for bridges between halos, e.g. Colberg (2007).) 
Around mean halos, the hlaments fan out and do not end 
within the scales (~ 800h“ *Mpc) we consider in 2 direc¬ 
tions. Increasing Ath > 0, these extra embedded structures 
decrease in size and eventually disappear, smoothly trans¬ 
forming into objects similar to those shown in Fig. 2. For 
Rs = 4Rr^Mpc this transition happens around Ath ~ 0.15, 
for Rs = 10h“*Mpc this happens at an only slightly larger 
Ath ~ 0.17. (As might be expected, this rich structure com¬ 
plicates using choice 1 to calculate A®, and is discussed in 
Appendix §B.) The scatter around this mean conhguration 
is signihcant as well, as discussed in §3.3 below. 

3.3 Mean profiles along axes and scatter 

Beyond shape, volume, extrapolated overdensity and mean 
Zel’dovich displacement of edges of each type of object, one 
can look at more details of the mean shear itself, for instance 
the shear values along any given axis. For the mean objects 
in Figs. 2,3, values along 3 axes for the 3 shear eigenval¬ 
ues (and the linearly extrapolated overdensity 5 = Ai) 
are shown in Fig. 6 . The dotted straight black line is the 
threshold Ath, crossing this means changing structure clas- 
sihciation (e.g. from a halo to a filament). The profiles for 


voids with A? corresponding to central shear choices 2 and 
3 (Ath = 0.762) are also shown, as these are explicitly cal¬ 
culable. Central shear choices 1, 2 and 3 correspond to line 
types dot-dashed, solid and dashed respectively. The differ¬ 
ent choices of shear at the origin give similar profiles for 
halos, but the filaments, sheets and voids show some large 
differences. The void profile, when calculable (i.e. choices 
2 and 3 at the origin), has mass compensation at a dis¬ 
tance away from the origin in the linearly extrapolated over 
density as expected**. Note that the mean shears take a 
very long time to approach their eventual value of zero even 
though they fall below the threshold Ath fairly quickly. 

One other scale in the system is the scatter, the vari¬ 
ance of the shear around the mean values (calculated in 
Appendix §A, Eq. A4).*^ The variance depends only on Rs, 
the window function, and the position relative to the con¬ 
strained origin, independent of the central value of the mean 
configuration and the threshold Ath- It is simplest to inter¬ 
pret along the 3 axes, where it starts at zero at the origin 
(where the shear constraint is imposed) and then increases 
to a constant ~ 0.2 at a distance ~ 27?s-|-l for one eigenvalue 
and at a larger distance, ~ 4.4i?a -|- 5.5, for the other two. 
(Calling the matrix in Eq. A4 Cij, where i,j are each pairs 
of indices, the scatter shown is \/{C~^)ii along i=ll, 22 
or 33 respectively.) The effect of this scatter on the classih- 
cation of a given region depends upon the configuration at 
the origin and Ath- We show in Fig. 7 the shear eigenvalue 
profiles along each of the 3 axes both for the hlament case 
seen in in Fig. 6 ), and for the filament configuration with 
Ath ~ 0 , but now also show the envelope of the variance. 

In both examples, far from the origin the mean value 
of shear approaches zero, while the scatter approaches its 
asymptotic value. At large distances, if Ath is high, the con- 
hgurations close to the mean do not scatter enough to change 
their classification from a void to anything else. As Ath drops, 
a one-sigma fluctuation will mix many different kinds of 
structures even at large distances from the origin. Thus, for 
instance,in Fig. 5 of the cross sections for low Ath, beyond 
about 20 h“*Mpc, the mean configurations are within one- 
sigma of all of the structures (e.g., the mean halo is within 
one-sigma of the mean filament, sheet and void). 

Thus the different ways of constructing the mean maps, 
via different constraints Co at the origin, provide, for the 
same Ath, similar objects in size and shape. The mean central 
objects increase in size and change shape as the threshold 
Ath is decreased, becoming very complicated as the threshold 
goes to zero. For very small Ath, the scatter strongly mixes 
the different structures even far from the origin. The scatter 
about the mean configurations is independent of the central 
value and approaches its asymptotic limit at larger distances 
from the origin as Rs increases. 


4 COLLAPSED HALO INTERPRETATION 

In general, any smoothing Rs, and any threshold Ath can be 
used to define halos, filaments, sheets and voids. One can 
also ask how objects in the smoothed map, such as halos, 

** We thank J. Peacock for suggesting we look for this. 

*^ We thank V. Desjacques for suggesting we include this. 
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Figure 5. The extreme behavior for Ath = 10~^ of mean regions around = Xi^cA (case 2) for a halo (top left) filament (top right), 
sheet (lower left) and void (lower right) for smoothing Rs = 4h~^Mpc. Colors are as in Fig. 2 and Fig. 3. The larger axes scales are 
used in order to capture the increased volume of non-void structures. The halos in filament edges and filaments in sheet edges visible for 
this Ath decrease as Ath increases, eventually disappearing for Ath ~ 0.15. 


are related to other ways of defining such objects (e.g. via 
gravitational collapse at later times). 


As halos are the most well studied collapsed objects, 
we investigated correspondences between our mean halos 
defined using Ai = ^i,CA and final collapsed halos such 
as are identified in numerical dark matter simulations. Ini¬ 
tial conditions for halos and their final properties have 
been extensively studied in Dalai, White, Bond ^ Shirokov 
(2008); Robertson, et al (2009); Ludlow and Porciani (2011); 
Elia, Ludlow & Porciani (2012); Despali, Tormen & Sheth 
(2013); Borzyszkowski, Ludlow &; Porciani (2014); Ludlow, 
Borzyszkowski & Porciani (2014). Final collapsed halos can 
be characterized by their mass and the measured linearly 
extrapolated overdensity of their initial regions. We checked 
whether we could find, for each Rs, a Ath such that our mean 
halo map overdensities matched ^extrap(Rs) from Despali, 


Tormen & Sheth (2013) and Vhaio Lag = M{Rs)/pb- We 
considered all three origin choices in §3.1.They gave simi¬ 
lar results: for smaller Rs, fixing {6) to the value matching 
Rs gave a Ath driven close to zero, or even below, and the 
Lagrangian volumes were much larger than that associated 
with M{Rs), and for larger Rs the opposite happened. For 
Rs = 4h“^Mpc, for cases 1,2,3 respectively we found that 
Ath = (0.37,0.48,0.43) and 1/ = (2.1,0.57,1.0)M(Rs)/p6. 

The quantities being compared are only approximate, as 
higher mass halos are also included, as higher peaks, with smooth¬ 
ing Rs , and there are different mass definitions for different quan¬ 
tities, with large conversion factors (White 2001). For the for¬ 
mer, the average mass of all halos with M > M{Rs) = M+(Rs) 
and their average 6+(Rs) = 1.686(1 -|- 0.2cr+(i?s) from inte¬ 
grating over the mass function, e.g.. Tinker et al (2008), gives 
Mj^{Rs)/M{Rs) > 1.2 until Rs increases to Rs 9h~^ Mpc. 
However S^{Rs)/S{Rs) ~ 1 for the smoothings Rs = 0.5 — 19.0 
which we considered. 
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Figure 6. Shear eigenvalue profiles Xi(r) (black, blue, red) and linearly extrapolated overdensity S (cyan) along x^y and ^ axes, for a halo 
(upper left panels), filament (upper right panels), sheet (lower left panels) and void (lower right panels). For the Ath as in Fig. 2,3 (and 
Ath = 0.762 for the void), profiles for our three cases of central shear are shown: (1) Aj = Xi^cA (medium weight, dot-dashed lines), (2) 
A° = k,CA (heavier solid lines) and (3) A9 = Aj^extr (lighter dashed lines). The straight black dotted line is the threshold Ath separating 
out the different structures (not visible for the void example at lower right). 
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Figure 7. Example of scatter around the mean profiles for a 
filament with two different central values, A? set using criterion 
1, and Ati, = 0.762 at top, using criterion 2 with Ath = 10“^ at 
bottom. Cross sections for their mean maps are shown in Figs. 2,5 
respectively. As a function of distance from the origin, as the 
smoothing Rs = 4:h~^Mpc for both, the scatter is the same for 
both examples. However, the effect on the classification of the 
regions differs because of the difference in the central values due 
to the difference in threshold Ath. The solid lines are the shear 
eigenvalue profiles Ai(r) (black, blue, red), and the filled regions 
show one sigma scatter along x, y and z axes. The top filament 
profile corresponds to the cross sections in Fig. 2 and the profile 
in Fig. 6, the lower to the filament cross section shown in Fig. 5. 
The dotted line again is the threshold Ath. The vertical solid red 
line is where the mean central filament transitions along that axis 
to a sheet. For low Ath, the scatter mixes different structures even 
at large distances from the origin. 


(This can be compared to the average mass for Rs and above 
at this scale, M+{Rs) ~ 1.5M{Rs).) Given the approxima¬ 
tions involved and the difficulty in taking into account the 
higher mass halos, this might be a close enough agreement. 
In summary, only for one smoothing scale Rs did it seem 
possible to find a Ath where the mean halo had its volume 
close to the Lagrangian volume M{Rs) / pb, for our 3 choices 
of mean object centers. 

Initial condition halos do not always give final ha¬ 
los, and not all halos at final times are tagged as ha¬ 
los in the initial conditions, for reasons which have been 
studied in detail in e.g. Dalai, White, Bond & Shirokov 
(2008); Robertson, et al (2009); Ludlow and Porciani (2011); 
Elia, Ludlow & Porciani (2012); Despali, Tormen & Sheth 
(2013); Borzyszkowski, Ludlow & Porciani (2014); Ludlow, 
Borzyszkowski & Porciani (2014). For instance, external 
tidal fields might destroy a initial local peak, or a peak might 
form at a later time due to a merger. 

If we take this choice of threshold as a preferred one 
for smoothing Rs = 4:h~^Mpc, because it gives a better 
connection to collapsed halos (of one fixed mass) at late 
times, then we also can consider the characteristic sizes for 
other objects for this smoothing, from Fig. 4. However, this 
lack of direct connection to collapsed halos (of one fixed 
mass) is not necessarily surprising. The mean objects we find 
for each kind of structure should be more characteristic of 
the smoothed map than, e.g., for halos, the collapsed objects 
which can form after full nonlinear evolution. 

For modified constraints on both Ath and S and taking 
A? = A i,CA, Papai and Sheth (2013) compared the ana¬ 
lytically predicted density anisotropy at large, i.e. close to 
linear, distances to that produced by a stack of numerically 
simulated halos, finding rough agreement. Properties of the 
mean objects for the origin constraints considered by Papai 
and Sheth (2013) are given in the Appendix §C. 


5 SUMMARY AND DISCUSSION 

A spatial distribution of Gaussian random fields can be clas¬ 
sified into halos, filaments, sheets and voids using a smooth¬ 
ing scale Rs and threshold Ath. These can be associated with 
linearly evolved initial conditions or with the initial condi¬ 
tions themselves for structure formation. Many smoothing 
scales and thresholds have been chosen for use in different 
contexts. The comparison of fractions of volume and mass 
in the different structures has shown some agreement be¬ 
tween analytic calculations of initial Gaussian random fields 
(specifically their linearly extrapolated values at the present 
time) and numerical simulations of the full nonlinear evolu¬ 
tion (e.g., Alonso, Eardley & Peacock (2014)). 

We explored further properties of cosmic structures in 
the Gaussian random field configurations, for the same clas¬ 
sification and parameters used by Alonso, Eardley & Pea¬ 
cock (2014). We used the properties of Gaussian random 
fields, combined with the power spectrum, to calculate prop¬ 
erties of mean halos, filaments, sheets and voids with three 
choices of central shears A^. These central shears are re¬ 
lated either to the mean shears for each type of object, 
or to the volume average of the extremal shear values for 
each object. For low Ath the mean structures are very com¬ 
plicated, while at larger Ath the mean structures centered 
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at the origin smoothly transition to the mean background 
through the sequence halo to filament to sheet to void. The 
mean halo has approximately similar extent in all three di¬ 
rections, the mean filament has one larger direction and the 
mean sheet is small along one direction and larger along 
the other two. We showed examples of mean structures at 
large and small Ath and discussed trends of their shapes 
with changes in threshold and in smoothing. One method of 
choosing the central shear was an implicit condition, which 
we could not always solve. When the central point is taken 
to be an extremal point (A® = Ai.extr), there is a natural 
value for the Zel’dovich displacement at the origin and thus 
one can also calculate the mean Zel’dovich displacement of 
the edge points of the mean objects, giving a final linear 
approximation to the size of each of these objects. For some 
Ath, the estimated Zel’dovich change in the mean filaments 
and sheets is small, while the mean filament and sheets are 
large relative to the smoothing scale (and sometimes the 
nonlinear scale), suggesting that these characteristic sizes 
might persist in the nonlinear maps in some form as well. 

We also characterized the scatter around the mean con- 
hgurations, which tends to a constant at a small multiple 
of Rs, the smoothing. The effect of scatter on the classifica¬ 
tion of structures around the mean background depends in 
part on Ath, we showed as an example the scatter for mean 
filaments for two values of Ath- 

There is not a clear identification of the mean halos 
with collapsed halos or collapsed halos above some mass. 
Matching the expected average linearly extrapolated over¬ 
density for a given smoothing Rs implies a threshold Ath, and 
leads to a Lagrangian volume (final estimated mass) which 
is too small (large) relative to M{Rs) for large (small) Rs, 
although there might be a sweet spot near Rs ~ 4h~^Mpc. 
These estimates were only approximate as different mass 
definitions are used in the quantities being compared. 

As these mean objects are an additional characteri¬ 
zation of initial Gaussian random field configurations and 
smoothed maps, it would be very interesting to compare 
them to mean objects in simulations with the same smooth¬ 
ing and threshold, both at initial and at late (fully evolved) 
times. This comparison is nontrivial. Although the division 
of halo, filament, sheet and void regions is direct in simu¬ 
lations, as are their volume and mass fractions, the separa¬ 
tion into specihc distinct objects, e.g. different hlaments, is 
not immediate. Splitting halo, filament, sheet and void re¬ 
gions into objects so that their stack gives our mean objects 
might be a natural way of defining such individual objects in 
the full map, but presumably additional assumptions will be 
required. The stacking might be most naturally associated 
with doing a volume average over extremal points as in our 
choice 3 of origin shear. 

This work is in a different direction than the extremely 
interesting research on finding the probability distribution of 
initial conditions which are consistent with a given final mea¬ 
sured nonlinear configuration (e.g. the line of study pursued 
in Jasche, Leclercq & Wandelt (2015); Leclercq, Jasche & 
Wandelt (2015a,b,c)) although both kinds of studies involve 
cosmic web classihcations of initial conditions and their in¬ 
ferred counterparts at late times. 

Several free parameters appeared in the construction, 
including the threshold Ath, the smoothing Rs and the cen¬ 
tral values A®. We chose a few different A®, however others 


are possible. Connections to other natural scales in the ini¬ 
tial conditions would be interesting to further pursue as well, 
such as those discussed in Shandarin (2009) and Pogosyan 
et al (2009). 

JDC thanks J. Peacock, A. Pontzen and J. Whitehouse 
for suggestions, and especially M. White for numerous dis¬ 
cussions. We thank N. Dalai, V. Desjacques, F. Leclercq, M. 
Neyrinck and the anonymous referee for helpful comments 
on the draft. 
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APPENDIX A: SHEAR CALCULATION 
DETAILS 

The mean shear {^ij (r)|Co) surrounding the origin with fixed 
boundary conditions {^ij I Co) = \%j can be calculated as 
described in Papai and Sheth (2013) (based in part upon 
Bardeen, Bond, Kaiser & Szalay (1986)). Their focus was 
on {(5(r)|Co), for only one of the boundary conditions Co 
which we consider, so we give more details for complete¬ 
ness. A more general angle-averaged probability distribu¬ 
tion P(^(ri),^(r 2 ),r) was found by Desjacques and Smith 
(2008). 

Defining^C(r)_^= (Cn (r), C 22 (r), ^33(r), ^i 2 (r), C23(r), 
^i 3 (r)), and if(0), Papai and Sheth (2013) use the Gaus- 
sianity of p(C(f)IC) to get 

(C(r)ia = ® (Al) 

which for conditions Co at the origin then give 

(C(r)|Co) = (C(r)a(C ® eraaCo) • (A2) 

The Zel’dovich displacement mean value, shown for central 
shear choice 3 in Fig. 3 can also be calculated analogously 
to that for ^(r) when the displacement at the origin is fixed 
to be zero (as it is for choice 3 of central shear): 

(4'f(r)|Co) = ('I'f(r)a(a®a“aaCo) . (A3) 

The variance of ^(r) — (^(r)|Co) is (Bardeen, Bond, Kaiser 
& Szalay 1986): 

(?® a - (a(r)a(a® a-^aao). (A4) 

which is zero at the origin (where there is a fixed constraint) 
and then increases, asymptotically reaching Q. 

These give the mean value of the shear, Zel’dovich dis¬ 
placement when Si'F = 0 at the origin, and scatter at posi¬ 
tion r when the shear at the origin is set by the boundary 


condition Co, which we explore in the body of this paper. 
The three central shear values A® which we consider are in 
§3.1. 

To evaluate Eq. A2, Eq. A3 and Eq. A4, we follow Papai 
and Sheth (2013) in choosing the shear at the origin to be 
diagonal ^(0) = ^ = ( 011 ^ 22 , 1 ^ 33 ), and substitute 

(O^WICo) = 

(6(0i(r)0i) - 3/2(02(r)02) - 3/2(0^(r)03))(0i|Co) 
+(-3/2(Oj(r)Oi) + 6(0i(r)02) - 3/2(Oj(r)03))(C22|Co) 
+(-3/2(Oj(r)Oi) - 3/2(0j(r)O2) J-6(0j(r)C33))(03iCo) . 

(A5) 

The expression for Jif can be found by replacing Oj (r) by 

'Pf(r). 

For the scatter, Eq. A4, 

((Oj(r) - (02(r)|Co))(04r) - (04r)|Co))) = 

(6(02(r)0i) - 3/2(02(r)02) - 3/2(0j(r)O3))(Oi|Oi(r)) 
+(-3/2(Oj(r)Oi) + 6(0j(r)02) - 3/2(Oj(r)03))(C22|04»’)) 
+(-3/2(Oi(r)Oi) - 3/2(0j(r)O2) + 6(Oi(r)C33))(C33iOi(»’)) 
+15((0i(i')O2)(O2OKr)) 

+ (0j(l')O3)(Cl3Ol(r)) + {irjOCAiCsiklO)) 

(A6) 

This is independent of the value at the origin, and the 
threshold Ath. 

The correlations (Oi(i’)Cmm) are 

+ '^3('>')(2fifmSjm + 2fj Cm^im (A7) 

-h , 

for mm = 11, 22,33, and 
1 z 

-(^?(r)Un) = i^mn -\-rmSin-\-fnSim) ■ 

(A8) 

The basis functions V['i(i') depend on the primor¬ 
dial matter power spectrum Ag{k) = CPL{k)/{2 'k‘^). 
We use a power spectrum for a AGDM universe with 
{Qrr,,^bC,h,as,n) = (0.31,0.022,0.68,0.83,0.96) (Planck 
Gollaboration 2013), calculated with the Eisenstein & Hu 
(1999) transfer function. For the 4'i(i') associated with ^ij 
(appendix of Desjacques (2008)), 

^ 1 W = A /o°° d\nkAl{k)WR{kfjAkr) 

’I'3(i') = A /o°° dlnkAj{k)WR{kf[-^j2{kr) - ^JAkr)] 

'&5(r) = ^ dlnkAj{k)WR{kf[7,joikr) 

+ SiMkr) + ^ji{kr)] . 

(A9) 

These are shown in Fig. Al for smoothing Rs = 4h ^Mpc, 
the case considered in most detail in this note. There are 
features at scales a few times Rs- For Jif’s correlations, the 
associated basis functions are (e.g.. White (2014)): 

’I' 2 (i') = d^nkAj{k)^^^J^{-lji{kr) - kj:^(kr)) 

’I'4(r) = dlnkAj{k) ^^^^^^ jaikr) . 

(AlO) 

These are used to calculate the mean displaced edges of the 
cosmological structures shown as solid black lines in Fig. 3. 
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Figure Al. The three kernel functions 'I's, ^5 (Eq. A9) used 
to calculate the mean shear properties, for a Gaussian window 
function with smoothing Rs = Mpc. The natural scale is 

about 10-20 h~^Mpc comoving, larger than the Rs in part due 
to the softness of the cutoff in the window function. 


APPENDIX B: SUBTLETIES IN FIXING 

A/ = Xi,cA 

Method 1 of assigning A° (in §3.1) constrains its value im¬ 
plicitly by requiring the average shear in a mean object of a 
certain kind, centered at the origin, to be equal to the mean 
value of shear for that kind of object in Eq. 11. In addition to 
not always being able to solve this constraint by systematic 
trial and error, there are also conceptual subtleties with this 
approach, although there is an intuitive appeal to requiring 
the average of the shear within the central mean object in 
our construction to be the analytically calculated average of 
the shear within that type of object, AqcA. One complica¬ 
tion was already noted for voids, as constructing a mean void 
volume is challenging when there is no boundary between a 
mean “void” configuration and the mean background, i.e. no 
mean void “size”. In practice, restricting to halos, filaments 
and sheets, we took objects centered on each kind of struc¬ 
ture and calculated the average shear within these objects 
for a given Ath- We then varied the central shear until the 
mean shear within each object was the one given by Eq. 11 
for that object. 

However, the mean maps often have more than just the 
central type of object and the void structure in them; all 
four structures can be present in a mean map centered on 


any one type of structure. So when one calculates the mean 
shear in the mean objects of a certain structure, one possible 
way is to include the contributions from all the different 
maps where that structure appears. That is all regions with 
that given structure, even if not part of the central object, 
should be included as well. If the average in the mean maps is 
chosen to include the contributions from these other regions, 
the calculation of mean shear becomes more complicated. 

We give an example by showing how adding these re¬ 
gions to calculate the mean shear can affect the calculation 
of mean shear for filaments. Generally, a mean halo at the 
origin also has filament regions around it. So in addition 
to the mean filament region around the origin, when the 
origin is a filament, the mean of filament shears including 
other regions would include the shear contribution from fil¬ 
ament regions around the mean central halo, calculated for 
the same threshold. The relative contributions of these fil¬ 
aments in the different maps, filaments around mean ha¬ 
los, around central mean filaments, (and sheets and voids in 
some cases), can be estimated if one breaks down the vol¬ 
ume fraction for each kind of object into contributions from 
each configuration: i.e. filaments when there is a mean halo 
at center, when there is a mean filament at center and when 
there is a mean sheet at center. This is 


^totfvf — Nhalo Vfilehalo “t“ NfuVfil Nsheet Vfilesheet 

“bNvoidsUilevoids ■ 


(Bl) 

Here the volume fractions fvh, fvf, fva, fvv of each struc¬ 
ture are available from Eq. 11, and Uiighaio is the filament 
volume associated with each mean halo, etc., found by con¬ 
structing the mean halo, etc. and measuring the volume. 
These volumes depends upon the central shears which one 
is trying to find. This is especially challenging if contribu¬ 
tions to the mean filamentary shear are large from any mean 
map besides that centered on the mean filament, and gets 
even more self-referential once filaments appear in the out¬ 
skirts of sheets or halos appear in the outskirts of filaments, 
i.e. the last two terms I4igaheet, Uflievoids 7 ^ 0 . 

For some Ath these additional terms NhaioUflighaio + 
NaheetUfiigsheet + Nvoids Ufligvoids Seem Small. Then the cal¬ 
culation ignoring these other terms, that is, averaging the 
filament shear only over the mean object centered on the 
origin which is a filament, seems justified even in this more 
complicated interpretation (and similarly for sheets and ha¬ 
los in their respective cases). For filaments these other terms 
can be estimated for size by the following method. For Ath 
large enough, the mean halo structure only appears when it 
is the central object, so A° can be found for the mean halo by 
imposing A° = AqcA for a mean halo centered at the origin. 
Also, again for Ath large enough, UfliGsheet = Ungvoids = 0. 
Then filaments only appear around the central mean halo 
and as the central mean filament. So we have: 


Vtotfvh — NhaloVhalo 

Vtotfvf = NhaloUfilghalo + NfllUfll . 


(B2) 


and the filamentary volume around the mean halo, Ehiehaio 
is determined by the mean halo, already in hand. As we also 
know fvh/fvf, we can constrain 

A'haloVfilGhalo _ VfilGhalo tpiqt 

AfllVai 

This can be used to see if the fraction of mean fil- 
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ament shear due to filaments around halos is large rela¬ 
tive to that due to a filament centered at the origin. For 
Rs = 4:h~^Mpc halos and Ath = 0.762, l/fiiehaio ~ 2Vhaio 
and fvf/fvh = 1/72, giving a fraction of filament volume 
due to filaments around halos of about 1/35, which we thus 
neglected in calculating the filament central value shown in 
Fig. 2. Repeating the argument for sheets, for this threshold 
and smoothing, the relatively large sheet to filament volume 
fraction fvs/fvf ~ 1/14, combined with the large Kheetefli, 
makes Kheetsfli too large to reasonably neglect in this inter¬ 
pretation.^^ In Fig. 2 we thus considered instead Ath = 1.25 
to impose At = \i,CA in a sheet configuration. (For the other 
scale studied in detail by Alonso, Eardley & Peacock (2014), 
Ra ~ lQh~^Mpc, filament contributions around halos and 
the filament to halo fraction (1/12) were both too large for 
their Ath = 0.261 to permit neglect of Etiighaio.) One could 
include these subleading contributions directly instead of go¬ 
ing to regimes where they can be neglected, as well, in some 
cases, when halos are not found in filament edges, and fila¬ 
ments are not found in sheet edges. It would be interesting 
to further explore this interpretation, although the issues 
with voids might make it difficult to push it very far. 

This mean halo/filament/sheet interpretation and ar¬ 
gument also gives a characteristic mean halo, filament and 
sheet volume and an expression for Afii in terms of the num¬ 
ber of halos, the volume fractions of filaments and halos and 
the mean volumes of filaments and halos. For the simplest 
case, for instance, Nm = ^ 44aio - Fkiehaio), and 

analogously for sheets. 


APPENDIX C: STACKED HALOS AND VOIDS 
WITH A DIFFERENT A? CONSTRAINT 

This method for calculating mean shear was used by Papai 
and Sheth (2013) to find the mean shear due to a given con¬ 
figuration at the origin. They took the origin to correspond 
to 1) the sum of all halos with mass above some smoothing 
scale and 2) all voids associated with a smoothing scale, at 
large distances and final times. To find the central values of 
shear, they imposed conditions both on the eigenvalues of 
the shear individually and on their sum (the density). We 
give here properties of the mean configuration maps for their 
boundary conditions at the origin. Note that we consider 
the mean smoothed values given their boundary conditions 
(that is we smooth the shear both at the origin and at the 
position it is being measured). 

They chose A® = Xi,Rs, where the latter was found 
by integrating Eq. 11 with either their halo or void con¬ 
straints but with an additional constraint in the measure, a 
step function in the density. For halos, Ath(o') = 0.41/cr and 
> 1-686/a, for voids, Ath(a) = 0 and J2i < —2.8/a. 
This is analogous to our choice of central shear 2, except for 
the additional density constraint. In addition, voids and ha¬ 
los have different thresholds Ath. 

We could also not find a A? (by trial and error) to meet the 
constraint At = Xi^cA by considering the volume centered on the 
sheet alone. More generally, configurations where any X^ qj^ < 0 
were often quite hard to find via trial and error. This happens, 
for Rs = Ah~^Mpc, at Ath 0.6 for filaments and Ath 1-2 
for sheets. 


For the comparison with simulations, as Rs increases, 
(j{Ra) decreases and Ath increases, so halos with a larger 
smoothing (higher mass) will also obey the constraints for 
a lower Rs (this was used in this note in §4). Thus for ha¬ 
los, Papai and Sheth (2013) interpret the resulting shear at 
the origin A? as due to the sum of halos with mass asso¬ 
ciated with Ra and above. They take simulated halos with 
masses corresponding to Ra and above in a numerical simu¬ 
lation, stack them, and compare measured and analytically 
predicted anisotropies in density at large distances, finding 
reasonable agreement. 

They consider masses 10^"^, 5 x 10^'*, Mq, from 

which we find Ra . In Table 1, we show halo and void config¬ 
urations they used, with the corresponding Rs, Ath, volume 
(~ M), shape, average z = 0 (Eulerian) overdensity, and A®, 
Ai. We use a Gaussian window function. 

Interestingly, for these constraints, the central value of 
shear and the average value in the halo or void volume are 
very close. The density cut is a stringent constraint on the 
shear and the average density in the regions is close to the 
minimum (maximum for voids) imposed density. The asso¬ 
ciated volumes for the halos correspond to smaller masses 
than the Lagrangian volumes of the halos associated with 
Ra. 


© 0000 RAS, MNRAS 000, 000-000 



Mean maps for cosmic web structures in cosmological initial conditions 


17 


M(Rs)/h-^MQ 

Rs/h ^Mpc 

Ath 

y Pb 
h-^MQ 

X, y, z max 

(5> 

Ai 

halo l.Oe+14 

4.2 

0.54 

4.70e+13 

5.52 6.67 4.36 

2.03 

1.35, 0.93, 0.65 
1.22 0.85 0.58 

halo 5.0e+14 

7.2 

0.81 

9.74e+13 

5.93 7.00 5.65 

1.79 

l.bd, 1.2(J, U.92 
1.53 1.14 0.86 

halo l.Oe+15 

9.0 

0.98 

1.03e+14 

5.52 6.36 6.22 

1.74 

1.84, 1.39, 1.09 
1.77 1.35 1.04 

void l.Oe+14 

4.2 

0.00 

1.20e+13 

4.06 3.13 2.59 

-2.88 

-0.77, -1.30, -1.84 
-0.74 -1.26 -1.77 

void 5.0e+14 

7.2 

0.00 

1.52e+13 

3.97 3.36 2.98 

-2.83 

-1.36, -l.yU, -2.42 
-1.33 -1.87 -2.38 

void l.Oe+15 

9.0 

0.00 

1.69e+13 

3.97 3.52 3.13 

-2.82 

-1.75, -2.28, -2.81 
-1.72 -2.25 -2.77 


Table Bl. Properties for halos (top) and voids (bottom) for objects using the smoothing scales Rs and thresholds Ath5<5 of Papai and 
Sheth (2013). These objects have both a Ath a^nd 5 constraint at the origin. For halos, Ath = 0.41/cr(i?s), <5 > 1.686/cr(i?s). For voids, 
Ath = 0, ^ < —2.S/o-{Rs). The mass of the mean halo volumes tends to be smaller than that associated with the smoothing scale. 
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